function[c,f,s]=pdefun(x,t,u,ux)  %建立偏微分方程函数
c=[1;1];

y=u(1)-u(2);

F=exp(5.73*y)-exp(-11.46*y);

s=F*[-1;1];

f=[0.024*ux(1);0.017*ux(2)];
